Determining role of Krein signature for three-dimensional Arnold tongues of 

oscillatory dynamos 
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Using a homotopic family of boundary eigenvalue problems for the mean-field cv 2 -dynamo with 
helical turbulence parameter a(r) — ao +^Aa(r) and homotopy parameter /3 € [0, 1], we show that 
the underlying network of diabolical points for Dirichlet (idealized, /3 = 0) boundary conditions 
substantially determines the choreography of eigenvalues and thus the character of the dynamo 
instability for Robin (physically realistic, f3 = 1) boundary conditions. In the (ao, /3, 7)— space 
the Arnold tongues of oscillatory solutions at f3 = 1 end up at the diabolical points for /3 = 0. 
In the vicinity of the diabolical points the space orientation of the 3D tongues, which are cones 
in first-order approximation, is determined by the Krein signature of the modes involved in the 
diabolical crossings at the apexes of the cones. The Krein space induced geometry of the resonance 
zones explains the subtleties in finding a-profiles leading to spectral exceptional points, which are 
important ingredients in recent theories of polarity reversals of the geomagnetic field. 

PACS numbers: 91.25.Cw, 91.25.Mf, 02.30.Tb, 02.40.Xx, 05.45.-a 
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I. INTRODUCTION 

Polarity reversals of the Earth's magnetic field have 
fascinated gcophysicists since their discovery by David 
and Brunhes [l| a century ago. While the last reversal 
occurred approximately 780000 years ago, the mean re- 
versal rate (averaged over the last few million years) is 
approximately 4 per Myr. At least two, but very likely 
three Q superchrons have been identified as "quiet" pe- 
riods of some tens of millions of years showing no reversal 
at all. 

The reality of reversals is quite complex and there is 
little hope to understand all their details within a simple 
model. Recent computer simulations of the geodynamo 
in general and of reversals in particular [H, 0, Q have pro- 
gressed much since the first fully coupled 3D simulations 
of a reversal by Glatzmaier and Roberts in 1995 Q . Most 
interestingly, polarity reversals were also observed in one 
Q of the recent liquid sodium dynamo experiments which 
have flourished during the last decade [3 Q . 

However, it is important to note that neither in sim- 
ulations nor in experiments it is possible to accommo- 
date all dimensionless parameters of the geodynamo [Ic| , 
and many of them are not even well known [Tij |. In an 
interesting attempt to bridge the gap of several orders 
of magnitude between realistic and numerically achiev- 
able parameters, Christensen and Aubert [l2| were able 
to identify remarkable scaling laws for some appropriate 
non-dimensional numbers. 

The use of appropriate simplified models [HI, [3 [TEI . 
[iH ] represents another attempt to understand better the 
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basic principle and the typical features of reversals. Most 
prominent among those features are the distinct asymme- 
try (with a slow decay and a fast recovery phase) [l7l ]. 
the clustering property of reversal events [181 ] , and the ap- 
pearance of several maxima (at multiples of 95 kyr) of 
the residence time distribution which has been explained 
in terms of a stochastic resonance phenonemon with the 
Milankovic cycle of the Earth's orbit excentricity (19LI20I]. 

One of the simplest reversal models which seems capa- 
ble to explain all those three reversal features in a consis- 
tent manner [2l|, [22| relies basically on the existence of an 
exceptional point in the spectrum of the non-self-adjoint 
dynamo operator, where two real eigenvalues coalesce 
and continue as a complex conjugated pair of eigenvalues. 
The importance of the specific interplay between oscilla- 
tory and non-oscillatory modes for the reversal mecha- 
nism had been early expressed by Yoshimura 12311 . Sar- 
son and Jones [2J|, and Gubbins and Gibbons [25|]. In 
the framework of a simple mean-field a 2 -dynamo with a 
spherically symmetric helical turbulence parameter a it 
was possible to identify reversals as noise-triggered relax- 
ation oscillations in the vicinity of an exceptional point 
[2H , I27I l28j . The key point is that the exceptional point 
is associated with a nearby local maximum of the growth 
rate situated at a slightly lower magnetic Reynolds num- 
ber. It is the negative slope of the growth rate curve be- 
tween this local maximum and the exceptional point that 
makes stationary dynamos vulnerable to noise. Then, 
the instantaneous eigenvalue is driven towards the ex- 
ceptional point and beyond into the oscillatory branch 
where the sign change of the dipole polarity happens. 
Therefore, the existence of an exceptional point is an 
essential ingredient for reversals, although non-linear dy- 
namics and the influence of noise must be invoked for a 
more detailed understanding of those events. 

From the spectral point of view, the reversal phe- 
nomenon of the geomagnetic field is strongly linked to 
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other fields of physics,hke van-der-Pol like oscillators 
[29| , geometric phases [3(| , 'PT-symmetric quantum me- 
chanics [sE [32j , PT-symmetric optical waveguides [33[ , 
microwave resonators [34j , and dissipation- induced insta- 
bilities us Hi nil . 

A particular problem of all those systems in which ex- 
ceptional points are involved is a strong sensitivity of the 
eigenvalues on boundary conditions (BCs). As for the 
geodynamo, the periodic occurrence of so-called super- 
chrons is usually attributed to the changing thermal BCs 
at the core-mantle boundary [ 2J, but the growth of the 
inner core may also play a role [281 ] by virtue of a spectral 
resonance phenomenon (38j . 

In this context it is worthwhile to note that impor- 
tant features of dynamos are most easily understand- 
able when treated with idealized (i.e. non-physical) 
boundary conditions. This was the case for explain- 
ing the famous eigenvalue symmetry between dipole and 
quadrupolc modes as it was done by Proctor in 1977 

Standing in this tradition, the present paper is devoted 
to a better understanding of the interplay of BCs, the 
spectral resonance phenomenon and oscillatory regimes 
in dynamos. 




FIG. 1: 1 = 0: (pink) spectral mesh |TT) for 7 = 0, /3 = 0; 
(dashed) eigenvalue parabolas (|12[) for 7 = 0, /3 = 1; 
(black) eigenvalue branches for /3 = 0.3, Aa(r) = cos(27rfcr), 
and (a) k = 1, 7 = 2.5, (b) k — 2, 7 = 3, with resonant 
overlaps near the locations of the diabolically crossed modes 
having (blue) the same and (green) different Krein signature. 



II. MATHEMATICAL SETTING 

The mean field MHD a 2 — dynamo [4l| in its kinematic 
regime is described by a linear induction equation for 
the magnetic field. For spherically symmetric a— profiles 
a(r) the vector of the magnetic field is decomposed into 
poloidal and toroidal components and expanded in spher- 
ical harmonics with degree I and order m. After addi- 
tional time separation, the induction equation reduces 
to a set of I— decoupled boundary eigenvalue problems 
[38| . which we write in a matrix form, convenient for the 
implementation of the perturbation theory [I^, IH| 



Lf 



0. 



l cHf + hd r f + l 2 f = 0, ilf 
The matrices in the differential expression L are 



(1) 



1, 



1 



h = d r l , 



1 



" ~ \-a(r) 1 

a(r) 



(2) 



and il := [21, 93] € C 4xS in the BCs consists of the blocks 



21 = 
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/3l + 1-/3 (3 
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(3) 



J g T f dr, where the overbar denotes complex conjuga- 
tion, and the boundary vector f is given as 



f := (f T (0),cU T (0),f T (l),d r f T (l)) e C 8 



(4) 



We assume that a(r) :— a® +jAa(r), where Aa(r) is 
a smooth real function with L Aa(r)dr = 0. For a fixed 
Aa(r) the differential expression L depends on the pa- 
rameters cxq and 7, while /3 interpolates between idealized 
((3 = 0) BCs, corresponding to an infinitely conducting 
exterior, and physically realistic ones [fi = 1) correspond- 
ing to a non-conducting exterior of the dynamo region 

The spectral problem {1} is not self adjoint in a Hilbcrt 
space, but in case of idealized BCs (j3 = 0) the funda- 
mental symmetry of the differential expression (38l . [lij 



L(A = 0) 



JL^J, 



1 

1 



(5) 



makes Lo selfadjoint in a Krein space (K, [., .]) [45[ with 
indefinite inner product [., .] = (J., .): 



[L f,g] = [f,L„g], f,ge/C. 



(6) 



For (3 7^ the operator Lo is not selfadjoint even in a 
Krein space. 



III. FROM DIABOLIC TO EXCEPTIONAL 
POINTS 



The vector-function f e H = L 2 (0, 1) ffi £2(0, 1) lives in 
the Hilbert space (H, (.,.)) with inner product (f, g) = 



In case of constant a— profiles a(r) = ao =const and 
j3 = the spectrum and the eigenvectors of the operator 
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matrix Lo are [381 ] 

K, = K ( a o) = —Pn + £aoVf>^ g 



e = ±, 



1 



fn G 



L 2 {0,1), neZ+ (7) 



with f n (r) being normalized Riccati-Bessel functions 
(2r) 1 /2j ;+i (^r) 

fn(r) = pf , 2 , s: , (fn'Jn) = Sn'n (8) 



and p n > the squares of Bessel function zeros 

Jl+ i (^/Pn~) = 0, < VPl < VP2 < ■ ■ 



(9) 



The eigenvectors f + , f„ € JC± C /C correspond to Krein 
space states of positive and negative signature e = ± 



(10) 



The spectral branches are real-valued linear functions 
of the parameter ckq with signature-defined slopes ±^/p~^ 
and form for all I = 0, 1, 2, • • • a mesh-like structure in the 
(a , 5iA)— plane. Spectral meshes for neighboring mode 
numbers I and / + 1 have only slightly different slopes 
of their branches and behave qualitatively similar under 
perturbations [38]. Therefore basic spectral structures 
for I = 1 , 2 dipole and quadrupole modes can be illus- 
trated by the simpler but unphysical / = monopolc 
modes which are given in terms of trigonometric func- 
tions. The (I = 0)— mesh built from p n = n 2 n 2 is de- 
picted as pink lines in Fig. [TJ 

The intersection of two branches A^, , A^ with n ^ n' 

occurs at points (ckq , A^) with 
™W — 



(11) 

and corresponds to double eigenvalues X^ — X n = X n , 
with two linearly independent eigenvectors and f*, , i.e. 
to so-called semi-simple eigenvalues (diabolical points, 
DPs) of algebraic and geometric multiplicity two [38], [43| . 

For the (I = 0)— mesh the diabolical crossings of the 
(en)th and (en+j)th modes with the same fixed \j\ G Z + 
are located on a parabolic curve [38| 



(12) 



A(a^) 



where ao — a — ir(2ns + j) and A^ 
ir 2 n(n + ej). Open circles in Fig. []Ja,b) indicate DPs on 
the parabolas \j\ = 2 and \j\ = 4. The Krein signatures 
[4(j| of the intersecting branches define the intersection 
index cr( u > = sS = sign (A'"-') in eqs. (jlljl . Branches of 
different signature 5 ^ e intersect for both signs of ao at 
< (green circles in Fig. [I]), whereas intersections at 
A^) > are induced by spectral branches of coinciding 
signatures: for e = S = + at ao > 0, and for e — S = — 
at ao < (blue circles in Fig. [T]) . 




FIG. 2: 1 = 0, Aa(r) = cos(47rr): (a) linear approximation 
of the 3D Arnold tongues and (b) their projection onto the 
(ao, 0)— plane indicating the influence of the intersection in- 
dex <t'"' on the inclination of the cones. 



For 7 = and constant a(r) = ao, the spectrum re- 
mains purely real on the full nomotopic family f3 € [0, 1] 
and passes smoothly in the (ao, 5RA)— plane from the 
spectral mesh at (3 — to non-intersecting branches of 
simple real eigenvalues for models with physically real- 
istic BCs at P = 1. For the monopole model I = the 
full spectral homotopy is described by the characteristic 
equation (1 — P)rj [cos (ry) — cos(a )] + 2/?Asin(?y) = 0, 
where 77(00, A) = ^Ja^ — 4A, which for physically realis- 
tic BCs (P = 1) leads to a spectrum consisting of the 
countably infinite set of parabolas (|12p labelled by the 



index j G Z + and depicted in Fig. Q] as dashed lines. 
The reason for the (P = 0)— DPs (fTTjl to be located on 
the (P = 1)— parabolas (fT2| is that the loci of the DPs 
are fixed points of the homotopy V/3 G [0, 1] — a phe- 
nomenon which indicates on their 'deep imprint' in the 
boundary eigenvalue problem (fTj). 

The eigenvalue branches with 3?A > 0, SA ^ (im- 
portant for the reversal mechanism [2(| HH) can be in- 
duced by deforming the constant a— profile into an in- 
homogeneous one, a(r) = ao + "/Aa(r), with simultane- 
ous variation of the BCs. This process is governed by 
a strong resonant correlation between the Fourier mode 
number of the inhomogeneous Aa(r) and the parabola 
index |j|. This is numerically demonstrated in Fig. [1] 
(black branches) for Aa(r) = cos(27rfcr) which highly 
selectively induces complex eigenvalue segments in the 
vicinity of DPs located on the parabola (|12p with index 
j = 2k. 

The underlying influence of 'hidden' DPs on real-to- 
complex transitions of the spectral branches can be made 
transparent by analyzing the perturbative unfolding of 
the DPs [H[ at the mesh-nodes (a^\x^) under varia- 
tion of the parameters ao, P, and 7. In first-order ap- 
proximation this gives for the (I = 0)— model 



X = X^-X^P+^-(a 



a, 



)± 



'D, 



(13) 




where an = n(en + 5n'), A^ — e5i: 2 n'n, and 

(en — Sn') (an — a^j 

i'ti \(el + 81)fA - (-l) n+n '(n + n')p7T 

(el - 51)fA - (-l) n - n '(n - n')pn "(14) 

with A := L Aa(r) cos[(en — Sn')irr]dr. 

For 7 = it holds D > 0, confirming that the eigen- 
values remain real under variation of the parameters an 
and P only. If, additionally, an = a§ , then one of the two 
simple eigenvalues (|13p remains fixed under first-order 
perturbations with respect to /?: A = )S V > in full accor- 
dance with the fixed point nature of the DP loci under 
the (3— homotopy. The sign of the first-order increment 
of the other eigenvalue A = A^ — 2A^/3 depends on 
the sign of A^ and, therefore, via (fTTj) directly on the 
Krein signature of the modes involved in the crossing 

In general, there exist parameter combinations yield- 
ing D < and thus creating complex eigenvalues. Eq. 
(|14p implies that in first-order approximation the domain 
of oscillatory solutions with SRA 7^ and SA 7^ in 
the (an , j3, 7)-space is bounded by the conical surfaces 

D = with apexes at the DPs (oq ,0,0), as shown in 
Fig. O Such domains, especially in case of r— periodic 
a— profiles, are in fact Arnold tongues corresponding to 
zones of parametric resonance [46[ in Mathieu-type equa- 
tions whose analysis in [47| was motivated just by Zel- 
dovich's studies on MHD dynamos. 

At the boundary D = the eigenvalues are twofold 
degenerate and non-derogatory, that is they have Jor- 
dan chains consisting of an eigenvector and an associated 
vector. Thus, DPs in the (an, P, 7)-space unfold into 3D 
conical surfaces consisting of exceptional points (EPs). 

The conical zones develop according to resonance selec- 
tion rules similar to those discovered in [38| for the case 
(3 = 0. For example, with Aa(r) = cos(2irkr), k £ Z, the 
constant A in (|T4"]) yields 



D := 



A = 



1/2, 2k = en-8n' 
0, 2k j= en - Sn' 



(15) 



so that in first-order approximation only DPs located on 
the (j = 2|£;|)-parabola (JT2J) show a DP-EP unfolding 
(in accordance with numerical results in Fig. [Tj). The 
cone apexes correspond to 2\k\ — 1 DPs with negative 
intersection index pd|) = — and countably infinite 
DPs with = +. The two groups are shown in Fig. [5] 
in green and blue, respectively. 

The real parts of the perturbed eigenvalues are given 

by 3?A = \ (u) (l-P) + a [ n V) (a - o^) I 2 and for fixed "o 
and increasing ft they are shifted (for both groups) away 
from the original DP positions toward the (5RA = 0)— axis 
(cf. the numerical results in Fig. [T]) — an effect which 
is similar to the self-tuning mechanism of field-reversals 
uncovered in 12611. 
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FIG. 3: / = 0: numerically calculated Arnold tongues for 
Aa(r) = cos(27rfcr), k = 2, and \ (v) < (green) or A M > 
(blue) and their approximations (dashed lines) (a) in the 
(00,7)— plane and (b) in the ((3, 7)— plane. 



Apart from this similarity, the eigenvalues of the two 
cone-groups show significant differences. The 2\k\ — 1 
cones of the first group have non-trivial intersection with 
the plane P = 0. In this (an, 7)— plane the zones of decay- 
ing oscillatory modes are 7 ^ —7 symmetric and defined 
by the inequality 



(an ± 2n(n - \k\ jf < 1 - 



1 - 



(16) 



where n = 1,2, . . . ,\k\. For k — 2 there are three primary 
Arnold tongues: 4chq < 7 2 and 16 (an ± 2ir) < 3j 2 . 
The cones of the second group meet the plane P = 

only at the apexes, having their skirts located in the 
sectors [p > 0,7 sign (an) > 0] and [P < 0,7 sign (an) < 

01 (cf. Fig. [2fb)). Therefore, in models with idealized 
BCs (P = 0) complex eigenvalues occur only in zones (fTB)) 
in the (ao, 7)-plane. 

The different oscillatory behavior induced by the two 
cone groups has its origin in the different Krein-signature 
defined inclination of the (D < 0)— cones with respect to 
the (P = 0)— plane. 

Passing from P = to a parallel (P 7^ 0)— plane, the 
(D < 0)— tongues (fTrj|) , corresponding to A^' < deform 
into cross-sections bounded by hyperbolic curves (black 
dashed lines in Fig.[3]Ja)) 

-4fc 2 (a ±27r(7i-|fc|)) 2 + n(2\k\-n)(~f ± 2Tr(n-\k\)P) 2 

> n(2\k\-n)4ir 2 p 2 k 2 , (17) 

with n = 1,2, ... , \k\. Since n < \k\, the lines 7 = ±2irnP 
and 7 = ±2ir(n — 2\k\)P, bounding the cross-sections of 

the 3D cones by the plane a^ = ±2n(n — \k\), always 
have the slopes of different sign. This allows decaying 
oscillatory modes for P = due to variation of 7 only. 

The (P 7^ 0)— cross-sections of the cones with (\W > 
0)— apexes have the form of ellipses (white dashed lines 
in Fig. Ha)) 



4fc 2 (a ± 27r(n+|fc|)) z 



+ n(2|fc|+n)( 7 ±27r(n+|fc|)/3) 2 
< n(2\k\+n)4Tr 2 p 2 k 2 , (18) 
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where n = 1,2,.... In the (/3 ^ 0)— plane the ellipses 
are located inside the stripe with boundaries 7 = («o i 
27r|fc|)/3 (pink lines in Fig. (Ha)), while the hyperbolas lie 

outside this stripe. Moreover, since in the plane o>q = 
±2n(n + \k\) the boundary lines 7 = ±2imj3 and 7 — 
±2ir(n + 2\k\)f3 have slopes of the same sign, the 7-axis 
does not belong to the instability domains, showing that 
for growing oscillatory modes the parameters (3 and 7 
have to be taken in a prescribed proportion, see Fig.[3Jb). 

The amplitude 7 of the inhomogeneous perturbation 
of the a-profile jAa(r) is limited both from below and 
from above in the vicinity of the DPs with a^ v > > 0. How- 
ever, numerical calculations indicate that this property 
can persist on the whole interval (3 € [0, 1], see Fig.[3][b), 
in agreement with the earlier findings of [48j ]. 



stantially determine the choreography of eigenvalues for 
(3 = 1 and, in particular, the loci of EPs which are im- 
portant to explain the reversals of the geomagnetic field. 
Although this has been exemplified for the unphysical 
monopole (I = 0) mode of a simplified spherically sym- 
metric a 2 dynamo model, the general idea is well gener- 
alizable to physical modes and to more realistic dynamo 
models. Work in this direction is in progress. 
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